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Abstract 

We present the results of relic density calculations for cold dark matter candidates coming from 
a model of dark energy and dark matter, which is described by an asymptotically free gauge group 
SU(2)z (QZD) with a coupling constant az ~ 1 at very low scale of Az ~ 10~^ eV while az ~ weak 
coupling at high energies. The dark matter candidates of QZD are two fermions in the form of 
weakly interacting massive particles. Our results show that for masses between 50 and 285 GeV, 
they can account for either a considerable fraction or the entire dark matter of the Universe. 
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1. INTRODUCTION 



It is almost universally accepted that the picture of the Universe made up of approxi- 
mately 4% baryonic matter, 23% dark matter and 73% dark energy represents a realistic 
cosmological model. However, it is astounding that almost 96% of the energy density 
of the Universe resides in some as-yet-unknown form. What is "dark matter"? What is 
"dark energy"? 

In Refs. [1, 2], a model of dark energy and dark matter was proposed in which a new 
unbroken gauge group SU(2)z - the shadow sector - grows strong at a scale ~ 10"^ eV. The 
gauge group SU(2)z was nicknamed Quantum Zophodynamics, or QZD, in Refs. [1, 2], 
where the subscript "Z" stands for the Greek word Zophos, meaning darkness. The model 
is described by an SU(2)z instanton-induced potential of an axion-like particle, Uz, which 
possesses two degenerate minima. The degeneracy is lifted by a mechanism described in 
Refs. [2, 3], yielding a false vacuum with energy density ~ (10"'^ eV)^ and a true vacuum 
with vanishing energy density. The present Universe is assumed to be trapped in the false 
vacuum [4], whose energy density mimics the cosmological constant. This is, in a nutshell, 
the dark energy model proposed in Ref. [2], which also computed various quantities of 
interest such as the tunneling rate to the true vacuum, etc. A Grand Unified Theory 
(GUT) involving the SM and SU(2)z was considered by Ref. [5] (The models presented in 
Refs. [2, 5] were later revisited by Refs. [6].). 

The particle content of the model includes two shadow fermions, i/'^f r) , with i = 1, 2, 
which transform as (1, 1, 0, 3) under SU(3)c<8)SU(2)L<8)U(l)y<8)SU(2)z, two messenger scalar 
fields (mediating between the QZD and SM matters; one of which is much heavier than 
the other [2]) with i = 1,2 transforming as ^1,2, = -1,3^, and one singlet complex 
scalar field (pz = (1/ 1/ 0, 1) whose imaginary part plays the role of the axion-like particle 
mentioned above. 

As discussed in Ref. [2], the masses of the SU(2)z triplet shadow fermions are found 
to be of the order of 100 - 200 GeV for the SU(2)z gauge coupling to grow strong at a 
scale ~ 10"^ eV, needed for the dark energy scenario. This coupling constant starts out at 
GUT-scale energy with a value comparable to that of the electroweak couplings, remains 
relatively flat until an energy comparable to the shadow fermion masses is reached, and 
then starts to grow after the shadow fermions drop out of the Renormalization Group 
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(RG) equations. At that dropout point, the SU(2)z gauge coupling becomes comparable to 
the weak SU(2)l coupling at the electroweak scale energy. These features have interesting 
consequences concerning the possibility of the shadow fermions being candidates for 
cold dark matter (CDM) in the form of weakly interacting massive particles (WIMP's) ^. 
The main reason is the fact that the annihilation cross sections for two shadow fermions 
into two SU(2)z "shadow gluons" are of the order of the weak cross sections, a typical 
requirement for WIMP's. The estimates that were made in Ref. [2] showed that it was 
possible for shadow fermions to be candidates for CDM with the right relic density. 

In this work, we would like to investigate this scenario in more details and by solving 
shadow fermions' evolution equations to determine the conditions under which they can 
be considered to be WIMP cold dark matter candidates. It will be seen that the mass 
range for the shadow fermions obtained by the requirement of having the right density 
fits in snugly with that used in the RG equations (i.e., the SU(2)z gauge coupling grows 
strong at a scale ~ lO""' eV). 

The outline of the paper is as follows. First, we go over the QZD model as far as the 
issue of dark matter is concerned. Then, we derive the evolution equations for shadow 
fermions and consequently solve them numerically, to obtain their relic density. Finally, 
the results of our relic density calculations will be presented and discussed, in comparison 
with the observational values. The shadow fermions relic density, when computed, would 
only depend on their masses. Therefore, the parameter space is simply two dimensional. 

2. THE SHADOW SECTOR AND ITS CANDIDATES FOR COLD DARK MATTER 

In this work, we only concentrate on the potential candidates for cold dark matter that 
QZD provides in the form of fermions. However, as discussed in Refs. [1, 2], the model 
offers a mechanism for leptogenesis through the decay of a messenger field, resulting in 
a net SM lepton surplus. 

For clarity, we list the particle content that is useful for our calculations, in particular 
the transformation of these particles under SU(3)c SU(2)l (8) U(1)y SU(2)z. 

• Two shadow fermions , with i = 1, 2, which transform as (1, 1, 0, 3). 

^ For a review on various features of CDM and WIMP, see, e.g., Refs. [7]. 
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• Two messenger scalar fields (p. , with i = 1,2, transforming as \1,2,Y^ = —l,3j. 
For relic density calculations, only the one with mass 0{< 1 TeV), i.e., <pf\ plays 
a role while the very heavy one with GUT-scale mass, i.e., is only useful for 
leptogenesis in this picture [8]. 

• One singlet complex scalar field (pz = (1/ 1/ 0, 1). The imaginary part Uz plays the 
role of the axion-like particle mentioned in section 1. The real part, Oz, was used as 
the inflaton in a model of "low-scale" inflationary universe [9]. 

We now briefly review the relevant aspects of the shadow sector that would be used 
in our relic density calculations for shadow fermions. 



2.1. The QZD Lagrangian 

The Lagrangian of Gsm ^ SU(2)z is given by [2] 

X = XsM + XL + Xvuk + Xcp -v{]^<f^^^f) - vljcpzf^ , (1) 
where Xsm is the SM Lagrangian and 

' i 

^^^^ = L L {s',Jl^?^'S + sipjL^pfVS) + Z ^'€^^^"5 + , (2b) 

m i 

£cp = J^G^[fv ■ (2c) 

In the above Lagrangians, G|^,^'s are the field-strength tensors of SU(2)z gauge bosons, 
the so-called shadow gluons, and the boldface typeset indicates the SU(2)z triplet mul- 
tiplicity. The sum over m is in fact over the number of SM families and the summation 
over i includes the number of shadow fermions. The coefficients gqi^^m, S<pimr arid K,- are 
complex. The covariant derivative in the Lagrangian can be written in the form 

D, = - z|r • W, - i^YB, - igzf ■ Af , 

where T 's are the generators of SU(2)z, which ought to be in adjoint representation when 
acting on shadow fermions, and A|f ^'s are the shadow gluon fields. The QZD Lagrangian 
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is invariant under a U(l)^ global symmetry, which yields an instanton-induced axion- 
like potential driving the present accelerating Universe. The transformations of QZD and 
SM particles under this U(l)^' global symmetry is given in detail in Ref. [2]. 

2.2. Masses and coupling constant 

The masses of shadow fermions come from the spontaneous breakdown of U(l)^\ 
Such a breakdown is made possible through the vacuum expectation value of (pz- There- 
fore, in the Yukawa coupling of shadow fermions with (pz, given in Eq. (2b), when (pz 
attains vacuum expectation value, ((pz^ = ^z, shadow fermions receive masses 

m^iz) = \Ki\vz, (3a) 

m^iz) = \K2\vz. (3b) 

The scalar messenger fields, on the other hand, are assumed to have zero vacuum 
expectation values to keep QZD symmetry unbroken. Their masses are non-trivially 
constrained by the evolution of QZD coupling, as explained in Ref. [2]. 

The QZD coupling constant, az = S^|/47i, is close to the SM couplings at high energies, 
while it increases to az ~ 1 at Az ~ 3 x 10"^ eV. The RG analysis of az, conducted in 
Ref. [2], studies the evolution of «z from Mqut to Az through a two-loop j6 function for 
possible masses of QZD particles. 

The RG analysis results indicate a direct correlation between the scale at which «z (E) 
starts increasing promptly and the mass of the lighter shadow fermion, nti. At energies 
prior to mi, «z (E) is mostly flat, but upon E ~ mi it begins to grow toward its value at Az, 
i.e., az (Az) ~ 1. 

Ref. [2] provides az (E) values for different conditions, i.e., masses, number of messen- 
ger fields, etc. However, a common thread among all analyses is that az does not change 
much from its value at Mgut until E ~ mi, being almost scale independent in that interval. 
At energies comparable to the masses of the shadow fermions, which themselves are of 
the order of he electroweak scale, «z is comparable to the electroweak SU(2)l gauge cou- 
pling. This will partially qualify QZD's shadow fermions as WIMP's and their candidacy 
for CDM, as already explained. 
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2.3. Shadow fermions as candidates for cold dark matter 

The two shadow fermions of QZD particle content meet the criteria for a WIMP, since 

• They interact very weakly with normal matter, i.e., through heavy scalar fields [2]. 

• They have cross sections of weak strength: masses in GeV and coupling constant in 
order of weak coupling [2]. 

• At least one is stable on cosmological scales: The lighter of the two shadow fermions 
is stable. The heavier one can decay into SM leptons and the lighter shadow fermion 
through the messenger scalar field (see Appendix B). However, if the shadow 
fermion masses are degenerate, both can be stable. Additionally, the shadow 
fermions can annihilate into shadow gluons or each other (if kinematically allowed). 

The messenger fields do not qualify as CDM candidates since they are unstable. The 
relic densities of shadow fermions can be obtained reliably by solving their evolution 
equations. Solving the evolution equations will reveal the applicable masses, which 
would give meaningful relic densities and put the model's candidates for dark matter 
into the test. 

3. EVOLUTION EQUATIONS FOR SHADOW FERMIONS 

The standard Boltzmann equation [10] describing the evolution of the number density 
n of a particle species i/^, is 



where H is the Hubble parameter, neq is the equilibrium density, v is the relative velocity 
in the annihilation process xl^xf> — > all, and (ov) denotes the thermal averaging of ov, with 
o being the total cross section of the annihilation reaction. The equilibrium density is 
given by 



where g is the species internal degrees of freedom and / (x, p) is the equilibrium distri- 
bution function. For particles that may play the role of CDM, the equilibrium number 




(4) 




(5) 
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density in the nonrelativistic approximation is 

3 

(mT\2 _,„ 

where T is the temperature, and m is the mass of the relic. The number density n 
satisfying Eq. (4) has two behaviors. In early times, n closely follows ngq but later when 
the temperature drops below m, the mass of the species, /leq starts to decrease exponentially 
until a"freeze-out" temperature is reached where the annihilation rate is not fast enough 
to maintain equilibrium. Below this temperature, n deviates substantially from and 
eventually gives the present day abundance of the species. Equation (4) can be solved 
numerically in relativistic (hot relic) or nonrelativistic (cold relic) regime. Ref . [11] showed 
that the validity of Eq. (4) and its solution breaks down if the relic particle is the lightest 
of a set of particles whose masses are near-degenerate and can contribute to the density 
of the relic through annihilation or decay processes, the so-called coannihilation case. 

For QZD's cold dark matter candidates, both shadow fermions can have present day 
abundances, if they have similar masses, which blocks the decay channel. For that reason, 
the evolution equations for the number densities of both species ought to be considered. 
The trivial reduction of shadow fermions occurs through their annihilations into QZD 
gauge bosons and the decay of the heavier one. Parallel to that, shadow fermions can 
annihilate into each other as well, which is analogous to the coannihilation case of Ref. [11]. 

To summarize, the reactions entering into Boltzmann equations for densities of shadow 
fermions are 

• Annihilation of shadow fermions into shadow gluons: i/'Pi/^P AS^^ A^^^ 

• Annihilation of a pair of one species into a pair of another: xjj^^^xjj^^^ ^ ^p^^'*-^^^\ 

• The decay of the heavier one into the lighter one and SM leptons: xpf'' — > U'lp^^K 

We assume negligible chemical potential for shadow fermions, which implies sym- 
metry among the number densities for particle and antiparticle of each species. To be 
inclusive, there can be a particle-antiparticle asymmetry in the shadow sector originating 
from the decay mechanism of messenger fields. The decay of messenger fields induces 
a particle-antiparticle asymmetry in SM leptons (see the decay of a messenger boson in 
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Fig. 11). The corresponding asymmetry in the shadow sector is expected to be as small as 
O (10"'') and therefore negligible to be considered in our relic density calculations. 

The evolution equations for number densities n\, of shadow fermions i/'^^^, 1/^2^^ are 
in the form 

+ 3Hni = -- {oxaO\a) (n\ - n^^^) - - (cTnt^ia) + 2 (^^21^^21) ^2 ' (6a) 

— + 3Hn2 = - - {o^AV^A) - n^^^^j - - <cr2iz^2i) "2 + 2 (^^12^^12) - T21 [n2 - n2,eq) , (6b) 

where r2i is the decay rate of the heavier shadow fermion, i.e., , Oij (with i, j = 1, 2, A) 
refers to the total annihilation cross section for the processes 

^(Z)^(Z) ^ ^(Z)^(Z)^ 

and Vij is the relative velocity of the annihilating particles for each reaction. Also, with a 
Maxwell-Boltzmann distribution function^, gq is given by 



(271) 

Tz 



grm]K2[f), (8) 



27l2' Viz- 

where Tz is the temperature of QZD matter, m, is the mass of the species and K2 is the 
modified Bessel function of second kind. The 1/2 factor on the right hand side of Eqs. (6) 
is to account for non-identical annihilating shadow fermions. 

Equations (6) can be written in a more convenient form by considering the number of 
particles in a comoving volume 

V, = f (9) 

which is the ratio of number density to entropy density, with the time derivative in the 
form 

dYi _ldni Hi ds 

If ~ lit ~ s^di' ^ ^ 

In the absence of entropy production (i.e., s = S/R^ with S = const.) 

ds ^ S 1 dR 



^ It has been shown that the use of correct statistics would only amount to less than 1% difference (see 
Ref. [12]). 
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which results in 



dYi drii 



The evolution equations, then, can be reformulated in the form 

dYi s 



dt 



= ^ [-{OiaVia) (Y\ - rteq) - <^12^12) Y\ + {o^^v^,) Y]\ , (13a) 



dYj s 



dt 2 



- (OiaVia) (YI - Yl^^) - <C72iy2l) Yl + {0^2Vl2) - -^21 (^2 - Y2,e^) 



(13b) 



s 

where y,;eq = ?^!,eq/s. Additionally, it is convenient to use the QZD plasma temperature Tz 
as independent variable, in place of time t. The relation between T (the photon tempera- 
ture) and Tz is easily found by the entropy conservation [1, 2]. The technique is essentially 
the same as that for finding the neutrino temperature using entropy conservation [10]. 
For example, at temperatures higher than the mass of the lighter messenger field (i.e., 
(fif'^) T > m^^, the QZD matter is in thermal equilibrium with the rest of the Universe, 
i.e., Tz = T. When T falls below the mass of the lighter messenger field, T < m^^, the 
QZD plasma conserves its own entropy separately and maintains its own temperature 
Tz T. The relation between T and Tz from there on can be found by entropy conser- 
vation anytime a particle decouples and transfers its entropy to the relativistic matter. 
At present, i.e., after e± decoupling, Tz = [(43/583)/(ll/18)]^^^ T. Ref. [1, 2] discusses the 
relation between Tz and T in more detail. Let us define = nti/Tz, we have 

dYj dYi dxj dYi ii^i dTz 



dt dxi dt dxi dt ' 



(14) 



where the time derivative of Tz satisfies 



'dTz\ ^ \ A ds 



(15) 



\ dt j SHsntidXi 
Considering all this, we can rewrite Eqs. (13) in their final forms 

^^''''''''''^ " J ^ ^''''''''^ " ^''''''''^ ' ^^^^^ 

dY2 X2 ds 



{02AV2a) [yI - Yl^^ + <C72lt^2l) " ((^^Vu) Y\ + -r2i {Y2 - Y2,,^ 



(16b) 



dx2 6Hdx2[ ^ ^'V 2 Aeq/ ^ - 2 X 1 g 

Equations (16) are first-order coupled differential equations in the form of Riccati 
equation, which ought to be solved numerically. The integration of Eqs. (16) from early 
Universe to present T° = 1.346 K (corresponding to photon temperature T = 2.725 K) 



yields today's number densities Y°. The present-day relic density of shadow fermion xp^^ 
in units of critical density pcrit is then 

a = — = - — - , (17) 

Pcrit Pcrit 

where Sq is the present-day entropy density of the shadow sector and pcrit = 3Hg/87iG. 
Finally, with Hq = 100/z km sec"^ Mpc"^ and Sq = 1271^7° ^45, Eq. (17) can be written in 
the from 

Qih^ = 0.5080 X 10^— , (18) 
GeV 

where h is the Hubble constant in units of 100 km sec~^ Mpc~^. Since ip^2^ decays, the 
relevant relic density is that of ip^^^- If ifii = m2, however, both shadow fermions can have 
present day abundances and only in such case, may we speak of two relic densities. 



4. COMPUTATIONAL METHOD 

Equations (16) include thermal averages {avYs, Hubble parameter H, and the derivative 
of entropy density ds/dxi, all of which need to be determined for numerical integration. 

The annihilation cross sections and the decay rate T21 can be calculated analytically. 
They are derived in Appendixes A and B and are given in closed forms, to leading order. 
The thermal averages (ovYs were then computed numerically using the compact integral 
form of Ref. [13]. In Appendix C, the relativistic thermal averages are provided in closed 
integral forms, expressed in terms of Xj. 

On the other hand, the Hubble parameter in a radiation-dominated Universe is given 

by 

H = ^^\nGp, (19) 

where G is the gravitational constant and p is the total energy density of the Universe, 
written as 

p = geff(T)^n (20) 

where geff 0^) is the effective number of relativistic degrees of freedom. Ref. [13] provides 
geii (T) values for two QCD phase transition temperatures Tqcd = 150 and 400 MeV. We 
made use of the geff i^) values corresponding to Tqcd = 150 MeV, which is a smoother 
function, as opposed to Tqcd = 400 MeV. It turns out that the solutions to Eqs. (16) 
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do not depend on the choice of Tqcd, mainly because the freeze-out temperatures for 
shadow fermions are always much higher than Tqqd, due to their large masses. As we 
already discussed, the relation between T and Tz can be easily determined by entropy 
conservation. As a result, the Hubble parameter in evolution equations was evaluated in 
terms of Tz and consequently x„ consistently. 

The entropy density s , in Eqs. (16), is mostly the entropy of the shadow sector. For 
temperatures T > m^p^, the QZD matter is in thermal equilibrium with normal matter and 
s is 

s = —g.sT', (21) 

where g*s = 459/4, and T = Tz- However, for most of the time T < and s is the entropy 
of the shadow sector, which is conserved separately, given by 



2n^ 
45 



Bosons Fermions 



(22) 



In both cases s is easily evaluated in terms of x„ providing values for ds/dxi of Eqs. (16). 

The numerical integration of the density evolution equations, Eqs. (16), was carried 
out using an implicit trapezoidal scheme^. We integrate from Xj = to Xj = nti/T^, 
where = 1.346 K is the present-day temperature of the QZD matter corresponding to 
To = 2.725 K, the photon temperature of the Universe today. 

Equations (16) were integrated for different sets of masses of shadow fermions varying 
between 30 and 300 GeV. The QZD coupling constant, az (£), values at energies Az < E < 
10^^ GeV are given for different sets of mi and m2 in Ref. [2]. Within the mass range we 
perform our relic density calculations, «z varies so slowly and continuously that it can be 
obtained for any set of and m2 by simple interpolation and extrapolation of the values 
provided in Ref. [2]. In this work, we have taken az dependence on mi, m2, and E into 
account in our relic density calculations. Nevertheless, it is worth mentioning that for 
a fixed m2 and at a given E, az does not vary much as mi changes. For example, from 



^ We implemented the idea of the backward differentiation formulas adapted to implicit trapezoidal 
scheme, presented in Ref. [14], for a system of Riccati equations. 
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Figs. (1-3) of Ref . [2], for m2 = 100 GeV and at E = 150 GeV one obtains 

az = 1.87500 x 10"^ or 4 = 3.51563 x 10"^ for mi = 1 GeV, (23) 
az = 1.87149 x 10"^ or al = 3.50247 x 10"^ for mi = 10 GeV, (24) 
az = 1.86567 x 10"^ or al = 3.48074 x 10"^ for mi = 50 GeV, (25) 

which demonstrate how az varies for 1 GeV < mi < 50 GeV . The az variation within 
such range (and similar mi ranges) is even less noticeable in relic density calculations, 
since we are dealing with «| in the annihilation cross sections. Our calculations showed 
that one could safely use an average a:| value over a wide range of mi values without any 
sensible loss of accuracy. For instance, an «| = 3.49961 x 10"^ for the above range works 
just fine. 

Ref. [2] carries out RG analysis of QZD's coupling constant considering a messenger 
field mass scale (mass of (pf'"' the lighter messenger field) niip^ = 300 GeV and higher, 
which points to when the QZD plasma decouples from the rest of the Universe. For our 
relic density calculations, we always chose m^^ > m2. It turns out that the relic density 
of shadow fermions does not depend on the choice of nifp^ > m2, as long as they are 
sufficiently apart^. That is mainly because the relic densities of shadow fermions (or more 
generally WIMP's) are mostly determined in their nonrelativistic epoch, i.e., for our case 
when Tz < m2. 

The decay of ip^2^ into a pair of SM leptons and ^pf'^ happens through a messenger field 
(see Fig. 11 of Appendix B). When the mass difference Am = m2 - mi is not very large, the 
decay rate for one of the possible decays can be given in an approximate form 

a^ tn^ 

^21 ~ TT? ^ (l - 8x + 8x3 _ ^4 _ ^2x^ ^\ ^ (26) 



28871 

where = g^J4n, F^^ is the decay width of the messenger field and x = m^lm\. As 
already said, we concentrate on the messenger field being sufficiently heavier than xp'"^^ 
where the "singularity" in the decay rate is not present, which can be seen from the 



* The thermal contact between the shadow and visible sectors may still be in effect through virtual exchange 
of a messenger boson for some temperatures below the mass of the lighter messenger field. With m^j 
being sufficiently larger than m2, the QZD plasma is ensured to have decoupled from the rest of the 
Universe before i/'j enters its nonrelativistic epoch, decoupling from an isolated QZD matter. 
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approximate from of T21, Eq. (26). We shall explain the interesting case of 1712 = in^^ when 
we present our results in the next section. It is worth mentioning, nevertheless, that such 
mass degeneracy poses no computational difficulty due to the presence of the messenger 
field's decay width F^^. 

On the other hand, is constrained for the model to predict the observed baryon 
asymmetry through an initial lepton asymmetry produced in the decay of messenger 
fields [2]. That requirement sets acp^ x 2.9 x 10"^^, which will consequently correspond 
to a long lifetime for \pf'^ (not less than 10^ sec). For that reason, the decay rate of xjjf'^ 
does not effectively enter the relic density calculations^, where the evolution equations are 
dominated by the annihilation processes. The remnant i/'2^''s (after the freeze-out) decay 
into SM leptons and j/'^^^'s anyway and we end up with no relic for i/^j^^ if the shadow 
fermion masses are not degenerate. 

5. RESULTS 

The relic density of shadow fermions depend on two parameters: their masses, mi, 
and m2. The masses affect the annihilation cross sections and consequently the dynamics 
of the evolution equations. Our relic density calculation results, therefore, are displayed 
either in terms of masses or mass difference. 

Suppose there were only one shadow fermion; in that case, the corresponding evolu- 
tion equation would be administered by shadow fermion's annihilation process and the 
expansion of the Universe. Since the annihilation cross section into shadow gluons and 
its thermal average (ov) are inversely proportional to the mass squared, a heavier shadow 
fermion would freeze out earlier than a lighter one, as it could not sustain a rate larger 
than the Hubble rate for as long. That would allow less time (at temperatures below the 
mass of the sole shadow fermion) for the Boltzmann factor to diminish the density, which 
would result in a higher relic density compared to a light shadow fermion's. This can be 
seen from the behavior shown by the dashed line in Figs. 1, 2, and 3, which describes the 
density of ij^f"^ or ^p'^^^ if they were the sole fermion in the QZD particle content. From 
those graphs, one sees that a heavy sole shadow fermion would have a higher relic than 

^ That means the decay of ipf''' is not determinant of the freeze-out temperatures. 
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FIG. 1: The relic density of 1/^2 versus xp^ s mass mi at fixed masses: solid lines, two shadow 
fermions at mi's (from left to right) = 50, 100, 150, 200, 250, 300 GeV; dashed line, one shadow 
fermion. Note that h = 0.732 in this figure and throughout this work. For the dashed line the 
horizontal axis is m. 

a light one. 

With two shadow fermions, however, there are two mechanisms governing the evolu- 
tion equations, besides the effect of the expanding Universe. There are those reactions, 
which exhaust the phase space from the species and those that populate it. The evolution 
of the number densities is determined by the competition of those mechanisms. The 
outcome of such a competition, on the other hand, depends on the masses of shadow 
fermions. 

Of those mechanisms, the decay of 1/^2^' plays no role in the early dynamics of evolution 
equations. Briefly, that is because the lifetime of which depends on mi, m2, and 
m^j, turns out either too long or too short to be a factor in the determination of freeze- 
out temperatures. For a well-separated set of m2, and m,pj, the lifetime of xp^^'' is within 
10^ sec < T2 < 10^'^ sec when mi m2, i.e, a nondegenerate case. That roughly corresponds 
to a temperature 1 keV < T < 1 eV, which is well after a typical freeze-out for i/^2^\ That 



means, the remainder of 1/^2^^ will decay into and SM leptons after the freeze-out, which 



(Z) 



leaves no present day abundance for 1/^2^^ The decay of an unstable shadow fermion at 
such low temperature into SM leptons can potentially disturb the cosmic microwave 
background (CMB). That, as we shall see, will place a bound on the mass of 1/^2^^ which 
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determines the density of \pf'^ at the time of its decay. With a mass degeneracy, i.e., 
mi = m2, of course ipf'^ is stable and decay is irrelevant. In that case, since the annihilation 
channel into another is also closed, we end up with two one-species cases: one for }pf^ 
and one for ip^^\ 

When m2 = m^p^, the decay width of the messenger field determines the lifetime of 
^p^^\ As discussed in Ref. [2], the requirement for the lightest messenger field to decouple 
before decaying yields F^j » m^^a^j, which is less than the expansion rate of the Universe, 
at T = m(pj . Since is of the order ~ 10"^'', a lifetime of 10"^^ sec < T2 < 10"^^ sec for 
is obtained. With such short lifetime, ^^^^ decays well prior to the decoupling of QZD 
matter, i.e., when QZD and the SM plasmas are in equilibrium. Effectively, that means 
we are down to the one-species case, regardless of the value of m2. 

Thus, the annihilation processes and their competition will mainly decide for the early 
dynamics of the evolution equations. At temperatures above the mass of the heavier 
shadow fermion ^^^^^ > both shadow fermions contribute to the population of another 
through the annihilation process i/^^i/^P — > ip^^^'ip^^\ As temperature decreases, the 
contribution of the lighter shadow fermion ip^^^ into the population of ip''^^ diminishes 
until it stops at an energy when it is no longer kinematically allowed. From there on, xjj''^^ 
will lose pairs monotonically due to its annihilations into shadow gluons and ip^^ pairs, 
while ip^^^ receives pairs from i/'j^^'s annihilation and at the same time loses pairs due to 
annihilation into shadow gluons. 

The annihilation into ip^-^^ provides an additional channel for ip''^^ to keep up with 
the expansion rate of the Universe and therefore delay the freeze-out. This reduces the 
density of xpf"^ prior to its freeze-out, compared to the one-species case, in two ways: (i) 
ipf'^ pairs are lost into ipf'^ pairs in addition to those lost into shadow gluons, (ii) the 
Boltzmann factor for temperatures Tz < m2 can act on density for a longer time. 

All this, though, depends on how apart i/'^ and 1/^2^^ are, masswise, at a fixed mi. 
Since the available phase space for ip^'^^ip''^^ — > ^p'f^Tp^^^ increases with the mass difference 
Am = m2 - mi, we expect i/'j^^'s density at freeze-out becoming small for an increasing Am 
due to a growing annihilation rate. On the other hand, a small mass difference reduces 
the phase space for the annihilation process and therefore increases the density. Knowing 
this is important in understanding the constraint on ^f'^'s density at the time of decay. 
Since the remaining 1/^2^' 's will decay anyway, there will be no relic for ip^^^ if mi m2. 
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which is reflective in Fig. 1, where s relic densities are displayed in solid lines for 
different m2's as mi varies. The relic density of ipf'^ falls down rapidly when the mass 
difference between the two shadow fermions is enough to allow the decay before our 
time and therefore to deplete the phase space from ip''^^ pairs. The maximum relic density, 
however, is always at nti = m2, where the annihilation cross section, Oij, and the decay 
rate T21 are vanishing and it is essentially the one-species case. 

The situation for \p^^^ is more complicated. The relic density of ipf'^ is shown through 
a solid line in Figs. 2, and 3 for different mi's as mi varies. For an extremely heavy 
ipf'\ i.e. mi = m2, 4^f'^'s relic density coincides with the one-species case, as expected. 
As Am deviates from zero ij^f'^ starts to dispense ij^'f^ pairs into the phase space (by 
annihilation earlier, and decay later) and thus Qi increases. Prior to freeze-out, this 
positive contribution comes from the pair annihilation of ip^^^ into i/'^^^, which will face a 
growing competition from ipf'^'s annihilation channel into shadow gluons, as mi declines. 
Since the annihilation cross section into shadow gluons grows for small masses, it will 
start to contend the rate of the extra ipf"^ pairs coming from 1/^2^' 's annihilation. For that 
reason, as mi decreases, the annihilation channel into shadow gluons depletes the phase 
space from ipf'^ pairs more effectively and therefore 4^f'^'s density before the freeze-out, 
which consequently diminishes its relic Qi. After the freeze-out, the remnant of ipf'^ will 
decay into ^pf'^ and lifts Qi, very much by a constant, except at small Am's where 1/^2^^ 's 
density is larger. 

For a nondegenerate mass case, i/'^^^'s relic density is what remains of shadow fermions. 
It is only at mi = m2 that the relic consists of both shadow fermions (equally so). To be 
inclusive of the degenerate case, the total relic density of shadow fermions Qj = Qi -I- Q2 
is presented in Fig. 4 against both masses and in Figs. 5, 6 against mi at fixed m2's, where 
the one-species case is also presented. The gray areas in Figs. 5, and 6 indicate the current 
bounds on the dark matter density from WMAP3 and all data sets [15]. 

It can be seen in Figs. 5, and 6 that the total relic density Qj increases as mi does, 
attaining a sharp maximum for the degenerate case, as if there were two "one-species" 
shadow fermions. On the other hand, Qj also increases with m2, which means for staying 
in the cosmologically allowed region a larger and larger mass difference would be needed. 
The two extremes are at m2 = 50 GeV, where the degenerate case is just making it to the 
allowed region, and m2 = 300 GeV, where a large mass difference is needed to stay 
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FIG. 2: The relic density, 0.ih^, of versus i/'^^''s mass nii at fixed ^p'"^^ masses: solid line, two 
shadow fermions; dashed line, one shadow fermion. For the dashed line the horizontal axis is m. 
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FIG. 3: The relic density, Qi, of versus ip^ s mass m-[ at fixed masses: solid line, two 
shadow fermions; dashed line, one shadow fermion. For the dashed line the horizontal axis is m. 
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FIG. 4: The three-dimensional depiction of the total relic density of shadow fermions as both 
masses vary. 

relevant. 

This can also be seen by looking at Fig. 8, in which the total relic density is displayed 
versus Am and the bounds are shown with two white dashed lines. We conclude that 
for m2 < 50 GeV, the total relic density is not enough to account for the total dark matter, 
even though the shadow fermions would still be relic particles taking on a fraction of the 
dark matter in the Universe. 

On the other hand, for m2 ^ 320 GeV, the total density of shadow fermions go beyond 
the upper bound and give unacceptable values even if we extend the mass difference to 
an extreme where mi = 1 GeV. That is shown in Fig. 7, where the total relic density at 
m2 = 318 GeV is only viable for a large mass difference of about 317 GeV and at m2 = 400 
GeV is no longer relevant. By going to such an extreme mass difference, we place a naive 
bound on the mass of the heavier shadow fermion, i.e., m2 ^ 320 GeV, above which the 
total relic density is no longer viable. 

There are, however, more restrictive bounds on m2 coming from the decay of ^p^^^ into 
SM leptons at low temperature, and its potential disturbance of the CMB of the Universe. 
We demand that 

1. The density of 1/^2^^ at the time of decay could not exceed that of the SM particles. 

2. The CMB density disturbance caused by the late decay of xp^^^ would not violate the 
CMB fluctuation, which has been observed to be at 10"^ level [15]. 
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FIG. 5: The total relic density of shadow fermions, Qrh^, versus i/^^^^'s mass mi at fixed jp'"^^ 
masses: solid line, two shadow fermions; dashed line, one shadow fermion. For the dashed line, 
the horizontal axis is the mass of the sole shadow fermion. The gray band represents the allowed 
density from WMAP3 and all data sets [15]. 
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FIG. 6: The total relic density of shadow fermions, Qj, versus \p-^ s mass mi at fixed 1/^2 masses: 
solid line, two shadow fermions; dashed line, one shadow fermion. For the dashed line, the 
horizontal axis is the mass of the sole shadow fermion. The gray band represents the allowed 
density from WMAP3 and all data sets [15]. 
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FIG. 7: The total relic density of shadow fermions versus i/'^ s mass nii at high xp^ masses. The 
gray band represents the allowed density from WMAP3 and all data sets [15]. 



Figure 9 illustrates these two conditions in graphs versus the mass of In Fig. 9 a, 
pi/pSM, i-e., the density of ipf'^ to the density of the SM matter - right before the decay - is 
plotted, which shows that the density of 1/^2^^ remains less than that of the SM particles for 
1712 ^ 285 GeV. The possible CMB density disturbance, dpyjpy, that the late decay of ipf'^ 
can create is shown in Fig. 9 b. The CMB density disturbance goes above the 10"^ order 
for i/'j^^s heavier than 245 GeV. The two above conditions, therefore, place a strong bound 
of 245 GeV on i/'j^^'s mass. 

As we discussed, the lifetime of ipf"^ could be very short if i/'j^^ and the messenger field 
were degenerate in mass. In that case, the total relic density of shadow fermions is simply 
that of the one-species case and it yields the right density for masses between 190 and 210 
GeV. 

For 50 GeV ^ ^2 ^ 245 GeV, the total relic density of shadow fermions can account for 
the amount of the dark matter in the Universe, depending on the mass difference. The 
total relic density lies within the observational bounds with small and even zero mass 
difference for light 1/^2^^ 's and with large mass differences when ip^^^ is heavy. 
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FIG. 8: The total relic density of shadow fermions versus their mass difference Am. The white 
dashed lines indicate the bounds from WMAP3 and all data sets [15]. 

6. SUMMARY 

We solved evolution equations for number densities of shadow fermions and obtained 
their total present-day density. The heavier shadow fermion turned out to be long lived if 
its mass differs from that of the messenger field. In that case, our results revealed an upper 
bound on the mass of the heavier shadow fermion, i.e., mi ~ 245 GeV, above which its 
late decay can potentially disturb the CMB density of the Universe beyond the measured 
fluctuation level of 10"^. 

For lighter shadow fermions, the total relic density can account for the entire dark 
matter of the Universe depending on the mass combination of shadow fermions. When 
the total density falls short of the observationally suggested density, it still, for most of 
masses, provides significant fraction of the dark matter of the Universe. 

Our results showed that if the heavier shadow fermion's mass is large, considerable 
mass differences would be needed to comply with experimental bounds. On the other 
hand, if the heavier shadow fermion's mass is small, little or even no mass differences 
suffice to give the right relic density. In that sense, degenerate and near-degenerate mass 
cases become relevant at low mass scales, but not for less than 50 GeV. 

A very short lifetime is expected for the heavier shadow fermion if its mass is the same 
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FIG. 9: Cosmological constraints on the mass of the heavier shadow fermion: a) The ratio of the 

(z) 

density of 1/^2 to the density of the SM matter, right before it starts to decay, versus the mass of 
ip^ , b) The amplitude of CMB density disturbance from the late decay of i/'2 versus the mass of 

/(Z) 

as that of the messenger field. In that case, the calculations reduce to a one-species case. 
Our results suggest that a sole shadow fermion must have a mass of about 190 - 210 GeV 
to account for the w^hole dark matter of the Universe. 

Last but not least, possible detections of the shadow fermion CDM candidates are 
briefly discussed in Ref. [2]. Needless to say, more work along this line is warranted for 
this model. 
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APPENDIX A: ANNIHILATION CROSS SECTIONS 

The pair annihilation of shadow fermions can yield either two shadow gluons or 
another pair of shadow fermions. The diagrams, to leading order, for both processes are 
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FIG. 10: Tree-level diagrams for: left, pair annihilation of shadow fermions into shadow gluons; 
right, pair annihilation of one type of shadow fermions into a pair of another type. 

displayed in Fig. 10, where the former process happens through three diagrams in t, u, 
and s channels and the latter in s channel. In those diagrams, p, p', k, k' , p„ p^, pj, p'. are 
momenta, I, I', n, n' and a, b are the QZD colors of shadow fermions and shadow gluons, s, 
s', Sj, s'^, Sj, s'. and A, A' are the spins of fermions and final polarizations of shadow gluons, 
and qt, qu, qs, q are momentum transfers. 

1. Annihilation into two shadow gluons 

We first compute the total annihilation cross section for a pair of shadow fermions into 
two shadow gluons denoted by three diagrams in Fig. 10. We carry out the computation 
for a fermion triplet with mass m, generically. The covariant amplitude A\ of the diagrams 
simply reads 



where, e.g., e„ (k) is the shadow gluon polarization four-vector, with A indicating its 
polarization state, qt = k - p, qu = k - p' , and qs = k + k' = p + p'. We are looking for an 
unpolarized cross section with the initial degrees of freedom averaged over and the final 




+ {qs-k') 



1^ T]ya + ik-qsX (p) , 



(Al) 
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ones summed over, which corresponds to the averaged squared amplitude 



A,A' s,s' a,b l,l',n 

We may even further compactify Ai in the form 



1 



(A2) 



(A3) 



with 



u-m 



m 



u^(p), (A4) 



where s, t, and u are the Mandelstam variables of the process. Therefore, lAif will have a 
compact form 

\Mf = ef (k) ef (k) el!' {k') e*^"' {k')%::;^'^% . (A5) 

The sums over the QZD colors of the squared amplitude, in Eq. (A2), result in five types 
of traces, namely 

^ (A6) 



a,b 

eacdebcdTr{f'r) = 12, (A7) 

a,b,c,d 

Y^ie,,Jr{rf''f'^) = 6, (A8) 

a,b,c 

Y^ieabJr(Pf'f') = -6, (A9) 

a,b,c 

Tr (ft^r^t") = 12 , (AlO) 

a,b 

knowing which yields 

colors 

+ 12 X 2Re (9Ji;atJ„) - 6 x 2Re (WJJls) + 6 x 2Re (aR*9Jit) ], (All) 

where the amplitudes SR^, ^Ru, 93is are colorless, having only the Lorentz degrees of free- 
dom. Evaluating lAif also includes summations over initial spins and final polarizations. 
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The sum over spins is simply the familiar y-matrix manipulation. On the other hand, 
sum over final polarizations involves terms like 

e"'^ (k) e*^"^ (k) and Yj ^ ) 

A A' 

To avoid closed loop diagrams containing ghost lines, we use the covariant form 



Y e^"' (^) e*"'' (^) = -'?^"' + 2 , (A12) 

A 

which preserves the gauge invariance and has the same effect as 

Y ^^'^ (^) (^) = -n^'' + ghost terms . 

A 

Considering all that, the spin averaged and polarization summed 9Jl- terms of Eq. (All), 
in terms of the Mandelstam variables of the process are 

'2 (w - m^) 8m"^ 



polarizations spins 



E ill"' 

polarizations spins 



2(t- m^) 4m2 Sm^ 



(A13a) 



u-nf" u-m" (u-m^) _ 
Tj I^^^I' = %[m'(2w-s)-m^-s'-t/(i/ + s)], 

polarizations spins 

L I L 2Re (aR;9Jl.) = [4m^ + - m^) + (. - m^)] , (A13d) 



(A13b) 
(A13c) 



polarizations spins 



polarizations spins 



polarizations spins 



s{t - m^) 



[m^ + {u -t)-{u + s)^] . 



(A13e) 
(A13f) 
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And finally, in terms of momenta, the unpolarized amplitude squared is given by 
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(p + pf\ P-k II 

The differential cross section in the center-of-mass (CM) frame, where p = (E, p) and 
p' = (E, -p), reads 
( do \ 



1 d cos 6 . 



CM 



nal\\\+v^cos^e ,\ l+u^cos^e 



3E2 V 



+2- 



1 - z;^ cos^ G 



(l-z;2cos2 



9\ UCOS0 
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, (A14) 



cos^ 6 

where v = |p| |e is the velocity of annihilating particles in the CM frame and az = gzj'^^- 



The total cross section then follows as 



na^ 1 
3E2 u 



2-v^ 



In 



1-v 



+ 



ln(l-v)---- 



In non-relativistic limit when E ^ m and v ^ 1, we obtain, neglecting O (v^), 

„2 



nr Z 



lOu _ 1 
2 



(A15) 



(A16) 



CM 3^2 V2z; 3 

The relativistic cross section in the lab frame (the rest frame of one of the annihilating 
particles) can be obtained as well. In terms of the velocity of the incoming particle in the 
lab frame v, it is 



2nal - (1 - v^) 



In 



+ 8i;2 + 4(2 - v^) Vl-z;2 - 8 
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2 - 2 Vl - u2 



(A17) 
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2. Annihilation of shadow fermions into each other 



The annihilation of shadow fermions into each other can occur through a shadow gluon, 
or the scalar field (pz- The smallness of the Yukawa coupling of (pz field, nonetheless, 
makes its channel rather negligible compared to the shadow gluon channel. For that 
reason and to leading order, the annihilation of a pair of ipf'^ with mass m, into a pair 
of ip^P with mass ntj is considered through the corresponding diagram of Fig. 10. The 
covariant amplitude Ai of the diagram reads 

M = (tl„„ (f "),„ «;/ (Pi) y,4 (p;) (p',) y""? (Pd , 

where, q'^ = s = + p'^j = (p, + p'i) ■ Once again, we are looking for an unpolarized cross 
section involving 



1 



V"" V"" i V"" V"" v 



Si,s'i Sj,s'. a,h 1,1' n,n' 

The gauge algebra calculations, which contain sums over QZD colors of the squared 
amplitude, result in a trace of the form 

Tr (f'T^) Tr (t^t") = 12. (A18) 

a,h 

After summing over QZD colors, we obtain 

I L \Mf = lxl2xm\ (A19) 

colors 

where the amplitude $R is colorless and only has Lorentz degrees of freedom. The Lorentz 
algebra including summations over initial and final spins yields 



spins 

^8' 



S2 



(A20) 



where s, t, and u are the Mandelstam variables of the process. The unpolarized squared 
amplitude is then given by 



■4 , 2 2 

(t-mj- mf^ + (ii -m^.-mfj +2 (m^. +mf^s , (A21) 
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and in terms of momenta 



1^1' = [(p/ • pO (p/ • + (pi • P') (Pi • pO + '^/P' • p^ + "^'Pi • Pi + ^^^y] ■ (^22) 

In the CM frame, where p, = (E, p) and = (E, -p), the differential cross section is 
I do \ 
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where = |p| |e is the velocity of the annihilating particles (i.e., ij^fipf) in the CM frame. 
The total cross section is then obtained as 
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The nonrelativistic limit of the total cross section, when Vi «: 1, and E — > m„ can be easily 
obtained, neglecting O {vf^, 
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In the lab frame, the relativistic total cross section can be also given as 
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where here is the velocity of the incoming particle (i.e., beam) in the lab frame. 
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FIG. 11: The decay of ip^ into SM leptons and ip^ through a scalar messenger field. 
APPENDIX B: THE HEAVIER SHADOW FERMION'S DECAY 

The decay of 1/^2^^ — > ll'xp^^^ is possible through the lighter messenger field cp^^^ (either 
real or virtual, depending on masses) and can yield any pair of SM leptons. Due to 
considerably small leptonic masses, when compared to shadow fermions', we carry out 
the decay rate calculation in the limit of massless SM leptons. In that sense, the decay rate 
for through cpf"^ with mass m,pj and a Yukawa coupling g(p^, representing any of g'^^„^, 
can be computed. The process, to leading order, occurs through the diagram of Fig. 11. 
The covariant amplitude Ai of the diagram is 

M = -gi,u^ (p), uf (p,), J I u- (pi), v^' (p\ , (Bl) 

where k = p2- p, T(p^ is the decay width of the messenger field and momenta pi, p2 refer 
to those of shadow fermions, while p, p' are the momenta of the SM leptons. Similar 
to previous cases, we are looking for an unpolarized decay rate with the initial degrees 
of freedom averaged over and the final ones summed over, which corresponds to the 
averaged squared amplitude 

^^ = \Yj\Yj 1^1'- ^^^^ 

Si,S2,r,r' jj 



^ 2 

There is not much of y-matrix algebra involved in computing \Ai\ , which easily gives 



„2 (pi • P') (P2 • P) 
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Finally, the decay rate, in the rest frame of ipf'\ can be found after the usual three-body 
decay kinematical considerations, which yields the integral form 



r ,(z) = do^—. 5 — (ml - ml) + (ml + ml) - 2a 

^2 72mlnJ^„^f ml - ml^ + m^Jl^ - cD^ 2 V V 2 ij 

(B4) 

where = and cD^ = (pi-pi) ■ The above decay rate behaves according to an 

dependence for W2 < mtp^ + m\, and an m^ dependence for m2 > mip^ + m\. 



m\ 



APPENDIX C: THERMAL AVERAGING 

The thermal averaging of ov (i.e., the annihilation cross section times the relative 
velocity) is discussed in Ref. [13], where a compact single integral for {ov) is provided. 
The authors of Ref. [13] explain that the thermal averaging of relativistic ov in the cosmic 
comoving frame and the lab frame are equivalent but they differ from the {ov) obtained 
in the CM frame. They stress that this difference is only significant in the relativistic limit. 
To stay relativistically covariant they introduce {ovmbi), for the cosmic comoving frame, 
where UmoI is defined in terms of the velocities of the two annihilating particles. The 
relation {ovmbi) = {ffVi^ah) {ffVcu) holds, in relativistic limit, anyway. To evaluate the 
thermal averages for our annihilation processes, we make use of the relativistically-valid 
single integral of Ref. [13], which is 

(ffUM0i) = T^jfr deVe(l + 2e)Ki(2xVrT7)c7z;iab, (CI) 
Kj (x) Jo ^ ^ 

where x = m/T {m is the mass of the annihilating particles), (x) is the modified Bessel 
function of order i and 



^lab = — T— r , (C3) 

1 + 2e 

with s being the usual Mandelstam variable for the annihilation process. The annihilation 
cross sections of shadow fermions are available analytically (see Appendix A). Therefore, 
the thermal averages of interest can be written with the help of Eq. (CI) in closed integral 
forms. The integrals then can be evaluated numerically for given masses. 
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For the annihilation of a pair of shadow fermions into shadow gluons ipf'^ipf'^ 
A^^^A^^\ the thermal average after simplification reads 

'e^ +4e + 2 



{OiAViA) 



Kl (xd 



I deKi [2xi VTTe) 

V/ 



(1 + e) 



3/2 



In 1 + 



1 + e, 



e^+3e + l 



In 1 



r~r~ \ 1 £3/2 5 - 

1 + e ) 6 1 + e 2 



(C4) 



(1 + 

where = nii/Tz- For the annihilation of one pair of shadow fermions into a pair of 
another, i/^^i/^P — > \p^^''xjj^^\ the corresponding thermal average is 



Iff: ]^ 



Jo ^ ' 1 + e \ mfl 



3mj Kl {xd 



+ e 



5e + 7e/6 



2 + ^+ 2 1 

6 1 + e 
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